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We first give a short review of the "local-current approximation" (LCA), derived from a 
general variation principle, which serves as a semiclassical description of strongly collec- 
tive excitations in finite fermion systems starting from their quantum-mechanical mean- 
field ground state. We illustrate it for the example of coupled translational and compres- 
sional dipole excitations in metal clusters. We then discuss collective electronic dipole 
excitations in Ceo molecules (Buckminster fuUerenes). We show that the coupling of the 
pure translational mode ("surface plasmon") with compressional volume modes in the 
semiclasscial LCA yields semi-quantitative agreement with microscopic time-dependent 
density functional (TDLDA) calculations, while both theories yield qualitative agree- 
ment with the recent experimental observation of a "volume plasmon" . 

1. Introduction 

Early in the history of nuclear physics, the (isovector) giant dipole resonance (GDR) 
provided one of the first manifestations of strongly collective excitations in finite 
fermion systems. Two classical models were suggested to describe the physics of 
the GDR: a) the model of Goldhaber and TellerjJJ in which protons and neutrons 
are both incompressible fluids undergoing a relative translational oscillation, and 
b) the model of Steinwedel and JenserP (previously also proposed by MigdaPl), 
in which protons and neutrons are both locally decompressed or compressed with 
opposite phases, such that the total nuclear density remains constant and a dipole 
oscillation results (cf. Fig.[l]below). Detailed analysis of experimental data revealed 
later tha t a suitable combination of both models was necessary to interpret these 
datajl^l^ so that the GDR could be classically best understood in terms of coupled 
translational and compressional dipole modes. These early classical models were 
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later refined by the so-called "fluid dynamics"^ and the "sum rule approach"^ based 
on the selfconsistent mean-field description of collective excitations in the random 
phase approximation (RPA) These semiclassical models were successfully used to 
describe collective excitations not only in nuclei, but also in metal clustersj^^l^^^ 
where a similar coupling between translational and compre ssional di pole modes has 
been shown to well describe the collective optical responseff^ESIIll 

In this paper, we review the "local current approximation" (LCA), which en- 
compasses both the fluid-dynamical and sum rule approaches and can be derived 
from a variational principle on the same footing as the RPA, and quote some of 
its results for metal clusters. We then apply the LCA to collective electronic exci- 
tations in Ceo molecules, for which recent experimenti^ have revealed a "volume 
plasmon" , a broad high-energy shoulder i n the p hoto-ionization cross section (other- 
wise dominated by the "surface plasmon"liSI12l) which again can be semiclassically 
understood as a compressional component of the collective dipole excitation. 



2. The local current approximation (LCA) 

The stationary Schrodinger equation H\v) = [T -\-V)\v) = E^\v) for a many-body 
system can be cast into the following well-known "equations of motion" ^ 

(oia[i?,a]|o> - fito,(oiaa|o) = o, (i) 

where the operators and are defined by 

ot|o>-l^.), ai^^> = |o), a|o) = o. (2) 

|0) is the ground state and hiUi, = E^ — Eq {i' > 0) are the excitation energies. As 
shown in Ref.^l^, Eq. ([1]) can be rederived by the following variational principle: 

SEslQySQ^O, EalQ] := ^m3[Q]/mi[Q] , (3) 
where the "moments" mi and ma - cf. ()15|) for their names - are defined by 

:=^(0|[Q,[^,Q]]|0), msiQ] -.= ^(01 [[H,Q],[[H,Q],H]]\0). (4) 

As long as Q is taken to be the most general (nonlocal) hermitean operator, the 
system p|4p is equivalent to the exact stationary Schrodinger equation. Successive 
orthogonalization of Qi, Q2t-- yields the exact excitation spectrum E^{Q^) = fUv^ 
(i^ = 1, 2, . . .). With Q,y (x Ol+0,y, which may be interpreted as a set of generalized 
coordinates (cf. Refs.lSII3|), we are brought back to 

In the selfconsistent microscopic mean- field approaches, one replaces |0) either 
by a Slater determinant (Hartree-Fock theory, HF) or by the Kohn-Sham (KS) 
ground state in terms of the local density p{r) (density functional theory, DFT). 
If the operator Q in ^ is replaced by a one-particle-one-hole (Iplh) operator, its 
variation ^ leads to the RPA equations. [The RPA ground state should in principle 
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contain 2p2h excitations, but due to a theorem by Thoulesi^^ one may use the HF 
ground state for computing mi[(5i/] and m^lQ,,] in (j4|).] In the framework of DFT 
(using the local density approximation, LDA, for the exchange-correlation energy), 
the RPA is often also referred to as the time-dependent LDA (TDLDA). 

In the sum rule and fluid dynamical approaches, one approximates the collective 
excitation energies hiUi, by the energy E^IQ^], defined as in ()3l4p but in terms of 
suitable local model operators (5i/(r), such as Qd = z to describe pure translations 
(Goldhaber- Teller model), the monopole operator Qq = for radial compressions 
("breathing mode"); Q2 = ^^120 for quadrupole oscillations, etc. 

The local current approximation (LCA) consists in the following assump- 
tions. One uses for |0) in ^ the uncorrelated HF or KS ground state, like in RPA 
or TDLDA, and takes the operator Q in (|3l4p to be a local function Q{y). For local 
and spin-less (e.g. Coulomb) two-body interactions V , one then obtains 

vn f fV^ 
mi[Q]=m^[n] = ^jn{T)-xi{v)p{v)d\, u(r) = - — VQ(r) , (5) 

where u(r) is a local displacement field which is proportional to the collective 
current, see ([T^ below. (This justifies the name of LCA.) Note that if V^Q(r) and 
hence V • u(r) is zero, then one has incompressible collective flow; otherwise the 
collective motion involves local compression of the fermi fluid. 

The moment TO3[(5] is a more complicated functional of u(r), of the spatial 
density p(r) — X^i^i ^nd the kinetic energy density t{t) — Y^^=i \^4>i{^)V 

(and possibly a current density according to the "current-DFT"l^^ in terms of the 
ground-state HF (or KS) wave functions (^i(r]. The variation 5Ez[Q]/ 5Q{y) — 
leads to fluid dynamical eigenvalue equations!^ 

= {hu:,f ^ p(r) u, (r) {3=x,y,z) (6) 

yielding the spectrum fiWi^ and eigenmodes Eq. ([6]) represents three coupled 

nonlinear fourth-order partial differential equations for Uj (r) , which in general are 
extremely hard to solve. Because of their dependence on the wave functions , one 
also speaks of "quantum fluid dynamics" which includes the effects of zero sound .1^ 
A practical way ("finite-basis LCA") to solve ^ approximately consists in 
expanding the operator Q{y) in a finite set of basis functions {Qp(r)}: 

M 

Q{r) = Y,cMr). (7) 

The variational principle (|3l4p then yields a set of M secular equations, whose 
characteristic equation is: 

det\Cpp, -{nuj,fBpp,\^0, (p,p',i/ = 1,2,...,M) (8) 

with 

Bpp' - {0\[Qp,[H,Qp,]]\0), 
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Cpp. = {Q\[[H,QpU[H,Q.p,],H]]\Q). (9) 

Solution of dHI yields the excitation energies hujy and the operators Qi/(r) creating 
the collective states The characteristic equation ([8]) is the basic equation of the 
LCA and looks similar to one version of the RPA equations. The basis set {(3p(r)} 
must be suitably chosen; here a good physical intuition for the considered collective 
motion is of great help. 

In a further approximation, one may replace the HF or KS ground-state density 

9 1 

by a selfconsistent semiclassical (extended) Thomas- Fermi (ETF) type smooth 
density, and use the ETF functional tetf[p\ (and possibly a corresponding func- 
tional for the current density) in computing the ingredients of ([9]). Although 
this approximation misses the quantum shell oscillations in the densities, it has 
proven to be sufficient for the evaluation of collective excitation spectra in many 
cases.mm 

The LCA is equivalent to the "generalized scaling model" ^ representing the 
system by a collective Hamiltonian 



Hcoii = 2 X! ^^PP' ^p'^p' + ^PP' ^p^p') (10) 

P:P' = 1 

that describes coupled harmonic oscillations with the velocity fields Vp(r,t): 

Vp(r,i) = dp(t)up(r) (11) 

and the local currents 

(r , t) = Pa (r , t) Va (r , t) , (12) 
which obey the continuity equation: 

^p„(r,t)+V-j„(r,t) = 0. (13) 

Here pa{r,t) are the "scaled" time-dependent densities (see Refs.^^^1^^^ for their 
definition) . 

Having solved either ([8]) or (fTO|) , one knows the eigenmodes of the system and 
can calculate its response to an external excitation operator Qe^t- To this purpose 
one defines a strength function: 

Sq^^Ae) = ^ |(HQe.t|o)p5(i? - nw,) , (14) 

whose energy- weighted moments mk{(^ext) become: 

/•oo 

rukiQext)-- £;'''5Q„,(i;)d£; = ^(nw,)'=|(i/|Qe.t|0)p. (15) 

The photo-absorption cross section a{Lu) in the long-wavelength limit becomes 

a{u;) = (4^a>/3c) Sd^piE^ha;) , (16) 
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where Qext — Qdip — ez is the electric dipole operator and to is the frequency of 
the external electric field. 

For applications in metal clusters, the following basis set of local operatorJl^l 
has successfully been used: 

Qp{r) = erPYLo{0). (p = 1, 2, . . . , M) (17) 

We consider here only dipole modes {L = 1); for p = 1 we then have the electric 
dipole operator Qi ~ Qdip — ez, while for p > 1 we obtain compressional dipole 
modes with V^Qp ^ 0. Hence, a finite set with Af > 2 of these operators will 
allow for the description of coupled translational and vib rational dipole modes. 
(Alternative basis sets involve spherical Bessel functions.!!^ 

It is illustrative to consider the mechanism of this coup ling for the classical limit 
of a metal cluster in the spherical jellium modehEHHISS] jjgj-e ionic density is 
taken to be a spherical uniform charge distribution with radius R: 

pi{r)^epa&{R-r), (18) 

where the bulk density po is chosen such that the integrated total ionic charge is 
opposite to that of the valence electrons. In the classical macroscopic limit [N —> 
cxd), we take the electron charge to be opposite to that of the jellium sphere: Pei^) — 
—pi{r) and neglect the kinetic and exchange-correlation contribution to the total 
energy, so that the system is entirely described by the dominating classical Coulomb 
forces. On the left side of Fig. [U we illustrate the pure translational vibration of 
the electrons against the ions, described by the single dipole operator Qdip — e.z. It 
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Fig. 1. Schematic picture of collective oscillations of electrons (dotted lines) against ions (solid 
lines) in the spherical jellium model for a metal cluster. Left: pure translation, yielding the Mie 
"surface plasmon" (cf. the Goldhaber-Teller mode in the nuclear GDR). Right: compressional 
dipole mode, leading to a "volume plasmon" (cf. the Steinwedel-Jensen-Migdal mode in nuclei). 



leads to the so-called "Mie plasmon" or "surface plasmon" with the energy 

Ez{Qd) = hiUM,e = ^h^e^N/mR\ (19) 

where N is the number of single-valenced (e.g. sodium) atoms. This is what corre- 
sponds to the Goldhaber-Teller model for the nuclear isovector GDR. On the right 
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side of Fig. [T] we sketch a compressional mode brought about by a suitable combi- 
nation of operators (|17p with p > 1, which corresponds to the Steinwedel- Jensen 
(Migdal) mode for the nuclear GDR. 

As shown in the equation ^ can be solved analytically in the above 

purely Coulombic classical limit, with the following interesting solution: 
When M > 1 modes are coupled, one of which has p — 1 (pure dipole mode) and 
all others have M — 1 different, but arbitrary real values p > 1, the spectrum always 
consists of one surface plasmon with frequency uj = uJMie and M — 1 degenerate 
volume plasmons with the frequency lo = u j„ni - 

This result could only be obtained in RefP^^ for specific examples. A mathemati- 
cal proof, for the general case of operators (jl7p with arbitrary L, is given in the 
appendix of this article. 

The frequency lo^oI is the bulk plasma frequency ujpi of the corresponding metal 



= LUpi = y/3eymrl , (20) 



rg being its Wigner-Seitz radius. This "volume plasmon" can also be brought about 
as a pure compressional mode by the radial (monopole) operator Qq = er^ (which 
is also used to describe the nuclear breathing mode): 



EsiQo) = hw^oi = ^3 h^e^N/mR^ = VihujKUe ■ (21) 

This mode can, however, not be excited by the external dipole operator because it 
corresponds to L = 0, and therefore does not couple to the electric dipole field. The 
surface plasmon with frequency ujMie then carries all of the dipole strength. 

In finite clusters with realistic smooth electronic densities Pe{r), and including 
the quantum-mechanical kinetic zero-point and exchange-correlation energies, the 
operators (|17p with p > 1 can couple to the electric dipole field, the degeneracy of 
the M — 1 volume plasmons is broken and their eigenmodes carry a finite amount 
of dipole strength.-J^ Experimentally, this manifests itself in a broad shoulder, or 
sometimes a small extra peak, of the dipole absorption cross section, located some- 
what below the energy of the bulk plasmon ([21]) . This volume plasmon had already 
been anticipated in the early TDLDA calculations of Ekardt.^23j 

In Fig. [2] we illustrate the above with some results for singly-ionized sodium 
clusters (left side: Najy, right side: Na^j^). The dots in the upper left part on each side 
give the experimental photo-absorption cross section.1231 The lower left parts give 
the LCA resultsp^ the vertical sticks correspond to the eigenvalues hwi, weighted 
by their percentage of the total dipole sum rule, and the solid curve is obtained by 
convoluting them with a Lorentzian to simulate continuum effects (as is customary 
also in discrete RPA and TDLDA calculations). On the lower right on each side we 
see the ionic structure,S3j as obtained in Car-Parrinello type molecular calculations 
in the cylindrically averaged pseudopotential scheme (CAPS). "^^^ We see that the 
cross section of Na^y, with a non-magic number 26 of valence electrons, clearly 
exhibits two main peaks which are due to the large average deformation of the ionic 
structure (as well-known also for the GDR in deformed nuclei). This deformation 
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is prolate, and hence the higher peak has roughly twice the strength of the lower 
peak. These two peaks here are both translational surface plasmons, split by the 
average deformation of the system. The volume plasmo n is hardly recognizable in 
this case. The cluster NaJ,^ has a magic numbeJm^ A^=40 of valence electrons 
and therefore a nearly spherical electron cloud (in spite of the non-sphericity of 
the ionic structure) . Here the response is dominated by one relatively sharp surface 
peak and exhibits a plateau in the high-energy shoulder which can be identified as 
a volume plasmon, in fair agreement with the LCA prediction. 



Na27+ 




2.0n 
1.51 

i.o: 

0.51 



2.0 2.5 3.0 3.5 4.0 



3.5 eV 




2.0 2.5 3.0 3.5 4.0 



Na41^ 




Fig. 2. Optic dipole response of NaJ-, (left part) and Najj^ (right part). Upper left: Experimental 
result.^'* Lower left: LCA result. -'^'^ Lower right: Ionic structure from CAPS calculation. ^^'■^^ 



3. Coupling of surface and volume modes in Ceo molecules 

We now turn to the optic response of Cgo molecules, the famous Buck- 
minster fuUerenes. They consist of 60 carbon atoms, each providing 
four valence electrons, so that iV=240 electrons can oscillate collec- 
tively against the ionic structure. Indeed, a giant resonance peak has 
been predicted^ and experimentally observecP^ around 20 - 22 eV, and in the 
framework of the jellium modeP^ been interpreted as a Mie surface plasmon. 

More recent experiments"'^^ on the optic response of Cgo molecules have focused 
at higher energies, and a broad shoulder was observed around 30 - 45 eV (see the 
circles in Fig. [3|) which was interpreted as a volume plasmon. A two-Lorentzian fit to 
the single photoionization cross section positioned the surface plasrnonat 22±0.1 
eV and the volume plasmon at 38 ± 2 eV. A TDLDA calculatioiff^I^ using a 
readjusted version of the jeUium modeP^ yielded, besides the main surface peak at 
22 eV, an extra peak around 42 eV (see the solid line in Fig. [3]), without however 
revealing the precise nature of the corresponding collective motion. In an ensuing 
debate, the interpretation of the volume peak as a co mpressional component of 
the collective motion was challenged,!^ and defended^Sl with reference to the well- 
known situation in metal clusters. 
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Fig. 3. Optic response of Cgo molecules. Circles: Experimental photoionization cross section. 
Solid lines: TDLDA calculation'^^''^* and dashed lines: present LCA calculation (see text for de- 
tails) of the optic dipole response. Both theoretical curves have been blue shifted by 5.5 eV and 
scaled to the height of the experimental surface peak. 



We shall presently corroborate this interpretation in the LCA. We use the same 
jellium model as in the TDLDA calculationp^ in which the ionic distribution of 
the Ceo molecule is replaced by a spherical shell with radius R and thickness A: 

pi{r)=po[Q{R2~r)~Q{r~Ri)], (22) 

with i?i = i?- A/2, i?2 = R+A/2 using the constants R = 0.354 nm and A = 0.153 
nm; po is chosen such that the integrated ionic charge is opposite to that of the 
valence electrons. 

To obtain a first rough estimate of the results to be expected, we use the same 
schematic classical model with Peir) — —pi{r) and neglecting kinetic and xc en- 
ergies, as discussed above. Furthermore we keep only the leading terms in A/R 
(although this is not really a small parameter). The energy E-^^Q) with the dipole 
operator Qd — ez then yields a Mie surface plasmon at 



EsiQd) = ficjAUe = yJh?Ne'^/imR^A = 2L4 eV , (23) 
while the monopole operator Qo = er^ yields a volume plasmon at 



Ez{Qa) = hw^oi = V?.nLOM^e = ^/ N / TuR^ A = 37 eV, (24) 

in surprisingly good agreement with the experimentally fitted peak positions. This 
is, however, a coincidence, since the volume plasmon obtained with the monopole 
operator cannot couple to the electric dipole field. Coupling AI dipole operators of 
the basis set ([T7|) with L = 1 and p = 1,2,..., M, we obtain the following results 
for the ingredients of the characteristic equation ([8]): 



fi2 Kf f?p+p'-2 

Bpp' = (PP + 2) , (25) 

TOD 



^..'=(- 3 \^PP' + -,i^-p-p'-2pp')y (26) 
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The equation ([8]) can then again be solved analytically at leading orders in A/i?; 
the calculation is similar to that given in the appendix. The resulting spectrum 
consists of one surface plasmon at hujsurf — 19.8 eV and M — 1 degenerate volume 
plasmons at hwyoi — 31.3 eV. Thus, the coupling of translational and compressional 
dipole modes is seen to shift both peaks towards lower energies, the volume peak 
by a larger amount than the surface peak. 

To be more realistic, we now use the quantum-mechanical KS ground-state 
density Pe(r) of Ref.^^, include the kinetic energy using the ETF functional tetfIp] 

1 

(up to 4th order, cf. Ref. ), and the xc energy in the LDA. The numerical solutions 
for the LCA spectrum then converge for AI > 8 coupled modes. This is shown in 
Fig. m where the sum-rule weighted LCA spectrum has been Lorentzian folded 
with a width of F = 5 eV to simulate continuum effects (which were included in 
the TDLDA calculation). Since the jellium model neglects the ionic structure, the 
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Fig. 4. Convergence of the sum-rule weighted dipole response of Ceo molecules in the LCA ap- 
proximation with respect to the number M of coupled modes with p,p' = 1, . . . , M. The curves 
for M = 8 (long-dashed) and M = 9 (solid) cannot be distinguished (see text for more details). 



collective peaks in both the TDLDA and the LCA calculations appear at too low 
energies, as is known to happen also for metal clusters.!^ Therefore, the TDLDA 
curve has been artificially blue shifted by 5.5 eV in Rem, we have done the same 
with our LCA curves. The converged result for ill = 8 corresponds to the dashed 
curve in Fig. [3] above. We note that the LCA volume peak now is located at 38 eV, 
exactly like that of the experimental two-Lorentzian fit. 

To account for fragmentation channels that were not measured, the results of the 
TDLDA calculations were rescalecP to fit the height of the experimental surface 
peak; we do the same here with our LDA results in Fig.[3l We see that both theories 
underestimate the dipole strength in the high-energy region including the volume 
peak, but give at least a correct qualitative result. The LCA result, through the 
explicit use of operators (|17p with p > 1, confirms the nature of the volume peak 
as due to compressional components of the collective electronic motion. 

Concerning the difference between the two theoretical results, for which the 
same jellium model was used, we note that in the TDLDA calculation P^ I "^^ ! both 
cr(n = 1) and 7r(rt 2) type valence electrons have explicitly been included in 
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the microscopic calculation of the linear respose; the a electrons were found to 
contribute most dominantly in the volume peak region. In our LCA calculation we 
have, however, used only one type of valence electrons, which might explain the 
difference of the results particularly in the region of the volume peak. 

4. Summary and Outlook 

We have briefly reviewed the local current approximation (LCA), a semiclassical 
approach which can be based on a general variational principle on the same footing 
as the RPA. In the LCA, the ground state of a finite fermion system is obtained 
in the selfconsistent mean-field approximation and the collective excitations are 
described by coupled local operators creating the local currents of the collective 
motion. This approach, which had earlier been successfully used for collective exci- 
tations in nuclei and metal clusters, has here been applied to the optical response 
of Ceo molecules. In recent experiments, a volume plasmon has been identified in 
the photoionization cross section at an energy ~ 38 ± 2 eV, whereas the dominating 
surface plasmon, already earlier known, was located at ~ 22 ± 0.1 eV. TDLDA 
calculations reproduce this result qualitatively after applying an ad hoc blue shift 
to compensate for the missing ionic structure in the spherical jellium model used. 
With the LCA we obtain very similar results as the TDLDA, using the same jellium 
model and applying the same blue shift, and further using a Lorentzian folding to 
simulate continuum effects. Since, by explicit construction of the coupled local ex- 
citation operators, the collective currents are known in the LCA, we can identify 
the nature of the volume plasmon in Cgo as due to compressional components of 
the collective electronic motion with respect to the ions. 

It would we worth while to corroborate our semiclassical interpretation by de- 
termining the transition currents both experimentally and theoretically in micro- 
scopic TDLDA calculations. Further improvement of the theoretical description 
should include the ionic structure of the Ceo molecules. While this might be too 
time consuming with purely miscroscopic methods, the LCA appears to be an ideal 
economic tool for this because the semiclassical nature of the collective electronic 
currents is little affected by the ionic structure, as known from the correponding 
results in metal clusters. 

Acknowledgements 

We are grateful to M. E. Madjet, R. A. Phaneuf, J. M. Rost and U. Saalmann 
for helpful comments and for providing us with their data used in Fig. [3l One of 
us (M. B.) acknowledges the hospitality of the Physics Department, University of 
Nevada at Reno (UNR), during a sabbatical visit, and thanks S. Kiimmel and P. 
G. Reinhard for clarifying discussions. 

Appendix: Coupling of surface and volume plasmons in jellium spheres 
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Here we give a pro of of the general solution of the characteristic equation ([8]) 
surmised in RefM^ for the classical limit with Pe{r) — ~pi{r) and neglecting ki- 
netic and xc energies. We use the set of operators Qp{r) = rPYLQ{9) with integer 
p = 1, 2, . . . , AI . The equation then becomes, for fixed angular momentum L, 



det Cn- 



/(L)det|AppH =0, (p,p' = l,2,...,M) (27) 



where /(£) is a factor independent of p,p' , and the matrix Appi is given by 

_ (2L + \)pp' + L{L + 1)(2L - p - p') - A (2L + \)[pp' + L{L + 1)] 

{I+P + P') 

and the eigenvalue A is the squared ratio of frequencies: 

A = {uj/ujyoi)'^ ■ 

We can rewrite the matrix (j28p in the following form: 

_ (2£ + l)[pp' + L(L + 1)](1 - A) - L{L + 1)(1 +p + pO 
il+P + p') 

Now take out a factor (2L + 1) from all rows to rewrite ((27)) as 

det \App, I = (2i + l)^^det \Dpp, \=0, 
where the matrix Dppi is given by 



D 



ppi 



{l-\){pp' + F)/{l+p + p')-G, 



and the constants independent of p, p' are defined as 

F = L(L+1), G = L(L + 1)/(2L + 1). 
We now consider two cases: 

a) L = {): Then F = G = 0, and the characteristic equation becomes 
Aet\{l-\)pp'/{l+p + p')\ = (1- A)^^det \pp'/{l + p + p')\ ^ Q . 



(28) 
(29) 
(30) 

(31) 
(32) 
(33) 

(34) 



Since for any p,p' > Q the determinant on the r.h.s. above is never zero, we get M 
degenerate solutions with eigenvalue A = 1, i.e. with the volume (or bulk) plasma 
frequency uj — uj^oi given in (PTjl . 

b) L > 0: In this case F and G are non-zero and the matrix Dppi has the form 



ppi 



/(l-A)(l + F)/3-G (1- A)(2 + F)/4- 
(l - A)(2+F)/4-G (l-A)(4 + F)/5- 
(1- A)(3 + F)/5-G (1- A)(6 + F)/6- 

V 



G 
G 
G 



(35) 



Notice that in each element, the first term contains the factor (1 — A) and the second 
term is the constant — G. We now replace the first row by the difference between the 
first and second rows, the second by the difference between the second and third. 
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and so on, until we reach the last row in which we do not change anything. The 
determinant, whose value is not altered by these manipulations, then becomes: 

(l-A)En {l-X)Ei2 ... {l-X)EiM 

{l-X)E2i (1-A)i;22 ... il-X)E2M 



det \D. 



pp' I 



{1-X)Emi-G {1-X)E^ 



M2 



G 



(1 — A) Emm — G 



(36) 



where the Eppi are linear expressions in the constant F. Only the in the last row. 



I 



I 
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the additive constant —G remains, while aU other elements now are proportional 
to (1 — A). The characteristic equation therefore becomes 

det\Dpp,\ = (l-A)(^^-i)det|i?pp,| =0, (37) 

where Eppi is the remaining matrix after removing the factor (1 — A) from the first 
M — 1 rows in ([36|) : its determinant is linear in A. We thus get M — 1 degenerate 
volume plasmons with eigenvalue A = 1, i.e. with lo ~ lj^oI again. The last eigenvalue 
is difficult to find in general. But when any one of the p values equals L, the last 
eigenvalue is found to be A = LG/F — LjilL + I), corresponding to the Mie 
plasmon with angular momentum L, i.e., llil = \/L/{2L + I) uiyoi (cf. Ref .^^^ . The 
generalization of this proof to a set of M arbitrary real values of p is straightforward. 
For the dipole case L = \ one gets the result stated in the paragraph above ([20]) . 
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